Two-particle decays of Be meson 
into charmonium states 



V.V.Kiselev 

Russian State Research Center "Institute for High Energy Physics" , 
Protvino, Moscow Region, 142280 Russia 

and 

O.N.Pakhomova, V.A.Saleev 
Samara State University, Samara, 443011 Russia 

Abstract 

The factorization of hard and soft contributions into the hadronic decays of Be 
meson at large recoils is explored in order to evaluate the decay rates into the S, P and 
D-wave charmonia associated with p and vr. The constraints of approach applicability 
and uncertainties of numerical estimates are discussed. The mode with the J/ip in 
the final state is evaluated taking into account the cascade radiative electromagnetic 
decays of excited P-wave states, that enlarges the branching ratio by 20-25%. 

1 Introduction 

After the first observation of meson by the CDF Collaboration at FNAL |]1| in the 
semileptonic mode with the J/ip particle in the final state, 

B+ ^ J/^l+ui, (1) 

one expects rather a significant, factor 20, increase of statistics with in the same mode after 
RunIL However, the uncertainties in the mass measurements are essential in the semileptonic 
channel, since the neutrino momentum is not detected directly. So, the two-particle decay 
mode 

^ J/V;7r+, (2) 

is the most prospective channel for such the measurements. Therefore, we need a qualitative 
theoretical modelling in order to predict the basic characteristics of (|]). 

The dynamics of B^ decays was studied in various theoretical approaches: the QCD sum 
rules and potential models operated with the exclusive decays and gave the estimates 
for both the branching ratios and the total lifetime summed over the exclusive modes, which 
is consistent with estimates of inclusive decays and lifetime in the framework of Operator 
Product Expansion (OPE) combined with the machinery of effective theory in the form 
of nonrelativistic QCD (NRQCD) |Q. However, the feature of two-particle decays, which 
are studied in the present paper 

B+ ^ cc[2=+iL,]7r+(p+), (3) 

is specified by rather a large recoil momentum of charmonium cc[2^+ilj], where s denotes the 
sum of quark and antiquark spins, L is the orbital quantum number running from to 2, 
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and J is the total spin of charmonium. Indeed, in the framework of potential models the ap- 
proximation of heavy-quarkonium wave-function overlapping for the calculation of hadronic 
form factors can be explored in the region, where those wave-functions are not exponen- 
tially small, i.e. if an amplitude under consideration is soft enough, and the nonperturbative 
modelling in the form of wave functions is reliable, while at large recoils the behaviour of 
vertices for the quarks entering the bound states is significantly modified due to hard gluon 
corrections. Then, the exponential decrease of quark-meson form factors is replaced by the 
power-like one at large recoils. In that case, one could factorize the hard and soft amplitudes 
0, which was recently explored in the description of two-particle decays of Be [0 as well as 
heavy-light mesons 0. 

While in the decays into the S and P-wave charmonia were considered, in the present 
paper we develop the same technique for a more accurate analysis including the D-wave 
states of cc, too. So, the approach is based on the fact that in the heavy quarkonium we can 
neglect the binding energy e in comparison with the heavy quark mass, since by the order 
of magnitude it is determined by the kinetic energy of heavy quark and antiquark inside 
the meson e ~ mq ■ f^, where v is the relative velocity of quarks, u -C 1, so that e <^ mq. 
Moreover, the region of soft part in the heavy-quarkonium wave-function is determined by 
the meson size, which is about r ~ with the quark momentum pq ~ mq ■ v. Then we 
can apply the nonrelativistic wave functions in the amplitudes, where the virtualities /x^ are 
less than {mq • v)^, while at large recoils in the decays the hard factors of amplitudes with 
virtualities greater than {mq ■ v)"^ should be described with account of hard gluon exchange 
between the constituents of heavy quarkonia. The perturbative QCD can be used for the 
hard amplitudes, if /i^ ^ ^qcd- We check these conditions of hard-soft factorization and 
estimate the uncertainties of numerical results by the variation of charmed quark mass in 
the limits constrained by the excitation energy of P and D-levels with respect to S-one. 

In addition, we factorize the matrix element of light quark current by the vacuum inser- 
tion. This approach and limits of its applicability are discussed in Then, we deal with 
the hard approximation for the four heavy-quark operator omitting possible renormalization 
effects. We take into account the perturbative corrections to the effective nonleptonic weak 
Lagrangian. 

The paper is organized as follows: In Section 2 we present the basic model assumptions 
and general formalism, while the analytic expressions for the widths of Be decays into the 
charmonium states with the pion are given in Section 3. We describe the input parameters 
and present the numerical results, too. Then we take into account the radiative electro- 
magnetic decays of P-levels in order to estimate the summed J/ipir yield in the B^ decays. 
In Conclusion we discuss the obtained results and their uncertainties. Bulky analytical 
expressions for the decays with p are placed in the Appendix. 

2 Hard-soft factorization in Be decays 

Neglecting the binding energy, in the soft amplitude we put the mass of Be meson mi equal 
to the sum of b-quark and c-quark masses mb + mc, and the mass of cc state m2 equal to 2mc- 
Then the heavy quark and antiquark inside the bound state move with the same 4-velocity, 
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so that in the accepted nonrelativistic approximation we can write down 

Pi Pb Pc 

^1 = — = — = — ' (4) 



Pl_ 


- H. 
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rric 


rric 



V2 = — = — = — , (5) 
m2 rric rric 

where pi^2 are the momenta of decaying and recoil heavy quarkonia, respectively, and pq are 
the momenta of quarks composing the heavy quarkonia. 

However, at large recoils specific for the decays of B'^ — > cc\^''+^L,j\Ti^{p^), the conditions 
of @) and (§) could be valid only if we take into account the hard gluon correction with a 
large momentum transfer 



where is the mass of vr or p. At the tree level as well as with soft gluon corrections the 
prescription of @j and @ would give zero matrix element, since a little smearing by soft 
gluons responsible for the formation of wave-function results in the exponential suppression of 
overlapping at large recoils. Numerically, the characteristic virtuality in the hard amplitude 
is equal to 1.0 — 1.2 GeV^ for the charmonium in the final state with m2 = 3.0 — 3.5 
GeV. We see that such virtualities are large enough for quite a reliable use of perturbative 
QCD. Moreover, a characteristic relative momentum of heavy quarks inside the bound states 
under consideration is about p ~ 0.6 — 0.7 GeV, and the ratio p'^/k'^ ~ 0.3 — 0.4 is quite 
a small parameter for the expansion. Thus, the kinematical conditions in the decays of 
— i> cc[2=+iLj]7r"'"(p"'") favor the application of hard-soft approximation with the accuracy 
about 30%. Another source of uncertainty is connected with neglecting the binding energy, 
and it is more essential. We will further test it numerically by the variation of charmed 
quark mass from 1.5 to 1.7 GeV. 

A general covariant formalism for calculating the production and decay rates of S-wave 
and P-wave heavy quarkonium in the nonrelativistic expansion was developed in [|l^. In 
this approach, the amplitude for the decay of bound state (be) possessing the momentum pi, 
total spin Ji, orbital momentum Li and summed spin Si into the bound-state (cc) possessing 
the momentum p2, total spin J2, orbital momentum L2 and summed spin S2 is given by the 
following: 



Aipi,p2) = E ^L,.5,.(qi)(^i^i.; ^1^1.1 JiJi.) X (7) 

dq2 



LizSiz 



(2j^)3 E '^L2zS2zM{L2L2z;S2S2z\J2J2z)M{pi,p2,qi,q2) 



where M{pi,p2, qi, ^2) is the hard amplitude of process with truncated fermion legs entering 
the initial and final mesons as described by the diagrams in Fig. |I|. Here "^LzSzi^) ^ire the 
nonrelativistic wave functions for the heavy quarkonia. 

With the accuracy up to second order in relative momenta qi and q2, the operators 
TsSziPyl) projecting the quark-antiquark pairs onto the bound states with fixed quantum 
numbers can be written in the form 

r5i5i.(pi,gi) = -f- — ( — fi - ii + mc)^i[ — i)i + ii -mb), (8) 

Arricmb rrii mi 



3 




Figure 1: The diagrams with the hard gluon exchange in decays of —>■ cc[2^+iLj]7r+(p+). 

where = 75 for the pseudoscalar initial state 5*1 = and ^1 = ^{Siz) for the vector one 
Si = 1, and 

r52S2,(P2,g2) = ^{^{—h + h-mc)M—h-h + mc), (9) 

^ 4m^ 1712 

where ^2 = 75 for the summed spin 5*2 = of recoil quarkonium, and ^2 = ^('S'22) for S2 = I- 
Here e{Siz^2z) denotes the polarization of vector-spin state. The color factor b"^^ j \f?) for the 
singlet states should be also introduced in the quark-meson vertices. 

Making use of projection operators in and we write down the hard amplitude 
M{pi,p2, qi, 52) in the following way: 

M(pi,p2,gi,g2) = Tr [rt(p2,g2)7^r(pi,gi)Oj , (10) 

where for the decay with the vr meson in the final state we have 

O, = OW + Otfl, (11) 

Ol^ = 7|^Hc/.a,^(^^^-^j 1^3(1 -75), (13) 

with the notations 

/ III I I I I 11 / III 

ii = — i>i + ^1 - ^3, f 2 = — p2 + h + h, f = — h pi + 01 - m- 

nil ^2 ^2 ^1 

The factor ai comes from the hard gluon corrections to the four-fermion effective weak 
Lagrangian. Expanding M{pi,p2,qi,q2) over small parameters qi/rrii and g2/"^2 at qi = 
q2 = 0, we get 

DM dM 

M{pi,P2,qi,q2) = M{pi,P2,0,0) + qia^ |qi2=0 + g2a^ Iqi2=0 + 

Cqia (Jq2a 

1 , 

+ 7;l2aq2f3^ ^ |qi,2=0 + --- (14) 



2 dq2adq'. 



2/3 



where the above terms correspond to various quantum numbers, first, with Li = L2 = 
for the transition between the S-wave levels, second and third, with Li = 1 and L2 = 0, 
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Li = and L2 = 1 for the transitions between S and P-wave states, fourth, with Li = and 
L2 = 2 for the transition between the S-wave initial quarkonium and the D-wave recoil meson. 
Thus, we can easily find that for the various orbital states, the soft factors in the amplitudes 
A{pi,p2) are expressed in terms of quarkonium radial wave-functions in the following way: 

f d'^q / 3 

^ ^2L^(q)gag/3 = ^—R"iO)e^pip,Lz), (15) 

where Eaip.Lz) is the polarization of vector particle, and £0/3 is the symmetric, traceless, 
and transverse rank-2 polarization for the spin-2 particle. The above wave-functions are 
represented as the products of radial and angular functions: = Yll^{9, (p) RL{r). 

For the ^Pi charmonium state we get 

E^"(P2,^2z)(li^2Z,00|l, J2Z) =£"(P2, J2Z). (16) 

For ^Pj{J = 0,1,2) states the summation over the quark spins and orbital momentum 
projections results in 



E e'^{p2,L2z){lL2zAS2z\J2,J2z)e'{S2z) = \ ^^e^^''^P2MP2, J2z), J2 = I, (17) 

S2Z,L2Z 



i_ 

V2n 

e'"'{p2,J2z), J2 = 2. 



For cc[ii:>2l one has 



E ^"''bs, L2Z){2L2Z, 00|2, J2Z) = £°^(P2, J2Z) 
Liz 



For cc[3Dj(j=i,2,3)] states we get 

E e''\v2.Uz)^(2L2z^S2z\J2.J2z)e\S2z) = 

S2Z,L2Z 

~\fi (l^°/3 ^p(P2, J2Z) - Vap ei3{p2, J2z) ~ V/Sp £^(^2, -^2z)) , -^2 = 1, 

{£aa{P2, J2z)erPpa' PIq""' + e/3a(p2, ^2z)£™pa' PW') , ^2 = 2, (19) 

£al3p{P2, J2z), J2 = 3, 

= -g-P + El^ , (20) 

and eai3p{p2, J2z) is the symmetric, traceless, and transverse spin-3 polarization tensor. The 
summation over the polarizations gives the following ordinary expressions W^'■ 



where 



1 

E ^a{.P2, J2z)e*p{p2, J2Z) = Vap, (21) 
J2Z=-1 
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^ 1 1 

Yl Sa/3(P2, J2z)ep^{p2, J2z) = 7^ [VapVpa + VaaVpp) - IzVocpVpa , (22) 
J2Z=-2 ^ ^ 
3 



J2Z =— 3 



77? 



"^T^aaP PrfP"ip ~l~ Par{Pj3(jP"tp ~l~ P arfP PpP'ya j (23) 
1 / 

yPaf3'P^r{P pa ^" 'Pafs'Pja'Pp'ri ~l~ PafiP^pParj 
'{'V a'yP f3r{P pa ~\~ Pa-yPlSaPpri ^" P a'y'P B pP arj 
ar{P pa ~l~ Pfi-yPaaPpr} ~l~ P fijPapPar) 

To the moment, we completely define the procedure for the calculation of matrix elements 
squared, that we perform by the use of analytic calculation system MATHEMATICA. 



3 Results 

Neglecting the tt meson mass, we obtain the following analytical formulae for the widths of 

B+ cc[2«+iL^]7r+(p+): 

QtT 777,2 (1 ~ ^) 

, 32 |i?2(0)P (1 + a:)3 0,0 
r 5, ^ r^.TT = — F^^}^ Ux' -2x + 3f, 25 

^ V ^1 - 256 ^R^(l + xf 
rr/? ^ V ^1 - 256 |i?^(0)P(l + ^)^ ..Ql 

r(..^^.2.)^^.M|)P|l±|^ (30) 

^ , 256 |i?''(0)P (1 + x)^ .0 ,0 

97r 1112 (1 ~ ^) 
o_ , 5120^|i?''(0)P (1 + x)^ 

6Tl 7712 (1 ~ ^) 

3^ X 8192^|i?''(0)P (1 + a;)7 , , 

37r 777,2 (1 — ^> 



where 



x = andF = a,2G^y,^J,^|i?,(0)ra?. 
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In numerical estimations we use the following set of parameters: 



\Rm\' 
\Rm\' 

mow 



1.27 GeV^, 
0.94 GeV^, 
0.08 GeV^ 
0.015 GeV^ 



rric 


= 1.5 GeV, 




= 0.14 GeV, 


rrib 


= 4.8 GeV, 


fn 


= 0.13 GeV, 


mi 


= 6.3 GeV, 


Vbc 


= 0.04, 


m2 


= 3.0 GeV, 


as 


= 0.33. 



The values of wave-function parameters are taken from [|T2[. Then we get the estimate of 
direct J/ip yield associated with the pion 



6.455 X lO'^^al GeV. 



(34) 



The decay widths into different charmonium states and the vr meson are presented in Table 
ID as the fractions of decay width for J /ipn^ , while the absolute values depending on 

the choice of ai are given in Table El. 



Table 1: The yields of charmonium states in hadronic two-particle decays of Be meson 
represented as the ratios. 
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he 
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2.7 
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XcO 




1.6 
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Xcl 


'Pi 


0.016 




51 


Xc2 


'P2 


1.4 




4.0 






5.4 




3.6 




'Di 


2.8 




3.8 




'D, 


0.053 




31 






2.4 




4.2 



An additional source of J/ip mesons is the two-particle decay of Be meson with p in 
the final state: B^ cc[^^+^Lj]p'^ . Calculating the decay widths T{B^ cc[^'^+^L,j]p'^) 
can be done in the same way as for the widths T{B^ — > cc[^^+'^Lj]7i~^). Indeed, we use the 
factorization of light meson current, so that in decays with p we incorporate the substitution 
f^p'^ — > mpfpe'^, where £3 is the p meson polarization. Taking into account the numerical 
values of fp = 0.22 GeV and = 0.77 GeV, we get the decay widths 1(5+ cc[^s+iLj]p+)^ 
which are presented in Tables ^ Cumbersome analytical expressions for the decay widths 
of —>■ cc[2s+iLj]p"'" are given in Appendix. 

For the sake of comparison with the results obtained in the ordinary technique with 
the wave-function overlapping, in Table ^ we show also the estimates, which were recently 
evaluated in [^3|, |15|. The analysis of two-particle hadronic Be decays with the charmed S- 
wave recoil- mesons in the final state was also done in |]l6l , where the estimates are similar with 
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Table 2: The widths of Be meson with the charmonium states in hadronic two-particle decays 
calculated in the hard-soft factorization in comparison with the results of wave-function 
overlapping technique 0, p!5[| . 
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O C 1 1 T 
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18 aj 


XcO 




11 aj 


Xcl 


'Pi 
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cc[25+iLj]7r), 10-15 GeV 



r(5e cc[«+iL,,]p), 10-15 GeV 



2.1 al m 
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0.57 al [0 

0.32 a? [0 
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0.28 al 



27 a? 
26 al 
60 

37 a? 
5.2 al 
36 a? 
124 al 
70 a? 
11 a? 
65 al 
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5.5 al 
6.0 af 
1.4 a? 
0.81 al 1^1 
0.33 a2 [|13| 
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0.58 af PI 



those of [|15], but slightly less numbers. Then, we see that our estimates with the charmed 
quark mass fixed at rric = 1.5 GeV are significantly greater than the values calculated 



in [|T^, 15|. Indeed, for the S-wave charmonium the increase due to the nonexponential 



fall off the wave-functions is about a factor 4 in the matrix element squared, while for 
the P-waves this factor reaches one order of magnitude. The reason for such the increase 
is quite transparent. So, following the coulomb analogy, we can expect that the velocity 
of heavy quark motion in the P-wave quarkonium is less than the velocity in the S-wave 
state (remember, f„ ~ a/n, where n is the principal quantum number). Then, the wave 
functions of P-wave levels have more soft behaviour than in the S-wave states, i.e. the 
relative momentum of heavy quarks is softer in the P-wave states, while the overlapping of 
quarkonia wave-functions at large recoils is displaced into the region of high virtualities, and 
the suppression is stronger for the P-wave levels. Thus, we would expect the above result on 
the significant enhancement of P-wave level yield in the hard-soft factorization approach. 

Another problem concerning for the factorization applicability is related to inherent un- 
certainties due to neglecting the binding energy in the charmonium states. Indeed, putting 
rric = mj/^/2 leads to zero binding energy for the S-wave states, while for the excitations 
under study it is about 500 MeV, i.e. it could be quite essential in the numerical estimates. 
We test this dependence on the binding energy by the variation of charmed quark mass in the 
range of 1.5 — 1.7 GeV. We find that this variation brings the significant uncertainties into 
the absolute values of widths under consideration, so that the variations are about 30-50% 
for the P-wave charmonia and greater than 100% for the D-waves. Nevertheless, we observe 
that the ratios of widths presented in Table |l| are quite stable under such the variation of 
charmed quark mass. The corresponding uncertainties are limited by 5 — 10%. This fact 
implies that the theoretical predictions for the ratios of two-particle widths in the B^ decays 
are quite reliable. Moreover, these ratios are close to the values obtained in [O, ffH]. 
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Since the J /ip meson is experimentally detected with a good efficiency in the decays of 3^ 
we compare the direct J /ip yield {B'^ J / ipix^ {p'^)) with the cascade one in decays with the 
radiative electromagnetic transitions of excited P-wave states into J/ip. The corresponding 
branching ratios of radiative decays are known experimentally So, 



Br(xco J/^l) = 0.007, Br(xci Jj^l) = 0.27, Br(xc2 ^ J/^^) = 0.14. 
Then we obtain 

r(^+ ^ XcO,ci,c27r+ ^ J/^7r+7) _ ^ 

and 

f(Bj^W^ — = P"' 

Thus, we see that the correction to the J/ip yield in the hadronic two-particle decays of 
Be due to the indirect mechanism with the P-wave charmonium is about 20-25%. In con- 
trast, analogous contribution in the semileptonic decays of Be is significantly less, and the 
corresponding fraction due to the P-wave charmonium is about 5% 0. 
Supposing ai = 1.1, we get 

Br(5+ ^ J/^7r+) + Br(5+ ^ J/^p+) ~ 2.8%, (37) 

and in the Be decays to J ftp with radiative transitions of Xco,ci,c2 the correction is equal to 

Br(fi+ ^ J/^/'7r+7) + Br(5+ ^ = 0.64%. (38) 



4 Conclusion 

In this paper we have considered the hadronic two-particle decays of Be meson with large 
recoils in the technique of hard-soft factorization for the matrix elements. This factorization 
is based on the physical separation of hard rescattering the constituents composing the heavy 
quarkonia (the scale of virtualities ~ 1 — 1.5 GeV^) from the soft binding of heavy quarks (the 
scale of virtualities ~ 0.3 — 0.45 GeV^). Hard factors can be calculated in the perturbative 
QCD, while the soft ones are expressed in terms of wave-functions and their derivatives at 
the origin for S, P and D-wave levels of heavy quarkonia. We have calculated the widths 
of Be — > cc[2^+iLj]7r"'"(p"^) decays for the charmonium in the final states (see Tables 
We have found that the results for the ratios of widths are quite stable under the variation 
of model parameters, while the absolute values have rather large uncertainties especially 
because of variation of charmed quark mass, that reflects the main systematic uncertainty 
due to the approximation of zero binding energy in the charmonium. We have compared our 
results with the potential model |1^, |T5[ operating with the wave-function overlapping. The 
relative momentum of charmed quarks inside the charmonium states (especially inside the 
excited P and D- waves, where the relative velocity of heavy quarks becomes less than for the 
S-levels) is rather small in comparison with the recoil momentum, so that the overlapping 
is displaced into the exponentially suppressed region, where the wave-function formalism is 
not reliable. In this region the hard gluon corrections replacing the exponential behaviour of 
quark- meson form factors by the power one, are significant. Thus, we expect valuable yields 
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of excited charmonium states in the two-particle decays of Be- This increase results in the 
additional 20-25% fraction of J/ipir inclusive branching ratio due to the contribution caused 
by the radiative electromagnetic transitions of P-wave charmonium states. 
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Appendix 

In this appendix we write down cumbersome formulae for the ratios of decay widths with 
the p meson with respect to the widths with the pion for various recoil charmonium sates. 
We introduce x — m2/mi, where mi,2,3 are the masses of Be, ccf^^+^Lj] and p, respectively. 
Then we get 
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A _ UA. (23^"^ - 17x3 - 5x2 + + 29) \ / m|\ 
\ 3^ (1 - x) (1 + x)2 (x2 - 2x + 3) j ^ l^^y ' 

^ y 2mJ (1 - x)6 (1 + x)2^ 
-158x^ + 183x'' - 4x3 - 119x2 + 130x + 40) j + O ' 

/ m| x^- 6x^ + 32x^- 48x3 - 11x2 + 78x- 10 \ 

r(5e ^ Xc27r) " ^=^1^+6;;^ (1-X)2(1+X)2 ) ^ ' 

TjBe ^ CC[^D2]P) _ „ 

T{Be ^ cc[^D,]7r) ~ ^'^^l'^^] 



2ml 3x6 _ 12x5 + 27x'' - 16x3 - 3x2 ^ jqOx + 49\ / 
3 +0 I ^ 



mj (1 + x)2 (x3 — 3x2 _|_ _|_ 5^2 ) I ^4 J ' 



r{Be ^ Cc[3Di]p) ^ 

r(Bc ^ cc[3Di]7r) '^^"P^i] 



1 + ^ I fl69x^°- 

2ml (1 - a;)^ (1 + (5^c3 - 22x2 + 41x + 8)2 ^ 
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1306x^ + 5334x^ - 13168x^ + 22638x'' - 25436x^ + 10336x'^+ 
14672x3 - 8271x2-1514x + 642) j + Q ^ 



4 I ' 



822x^ + UOlx'^ - 124x3 - 873a;2 + 1230a; + 122) ) + O | ^ 



r(Be ^ cc[3D3]7r) " \^ + 12m? (1 - a;)^ (1 + x)^ + ^^"^ 



96a;3 - 15x^ + 152x - 22) ) + O ( ^ 



m 



The factors B for various charmonium states tend to the ratio of leptonic constants of p and 
TT mesons, if we neglect the p meson mass, while in general case we obtain 



x^ (1 - xf r 2^ 



(i+x) (i-x)2-3 



m 



1 , 



5.. = I — —Jii-x)Hi+xr-'-^ii+x^) + ^,, 

^'^(l+.)((l-.)^-"^3^ V 



e2 



B 



_n {i-xf 1 



„ _fl {l-xf{l + x) I 

-°cc[lD2] f2 / n\4 / — 2 T 
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^ccpD2] — ^ccl^D-i] — -^ccpDi]- 

Finally, we note that the most significant numerical correction due to the nonzero p meson 
mass comes from the phase space in the coefficients B. 
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